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Context. The understanding of fossil fields origin, topology and stability is one of the comer stones of the stellar magnetism theory. 
On one hand, since they survive over secular time-scales, they may modify the structure and the evolution of their host stars. On 
the other hand, they must have a complex stable structure since it has been demonstrated by Tayler and collaborators that simplest 
purely poloidal or toroidal fields are unstable on dynamical time- scales. In this context, the only stable configuration which has been 
found today is the one resulting of a numerical simulation by Braithwaite and collaborators who have studied the evolution of an 
initial stochastic magnetic field, which is found to relax on a mixed stable configuration (poloidal and toroidal) that seems to be in 
equilibrium and then diff'uses. 

Aims. In this work, we thus go on the track of such type of equilibrium field in a semi- analytical way. 

Methods. In this first article, we study the barotropic magnetohydro static equilibrium states; the problem reduces to a Grad-Shafranov- 
like equation with arbitrary functions. Those latters are constrained by deriving the lowest-energy equilibrium states for given invari- 
ants of the considered axisymmetric problem and in particular for a given helicity which is known to be one of the main actor of 
such problems. Then, we obtain the generalization of the force-free Taylor's relaxation states obtained in laboratory experiments (in 
spheromaks) that become non force-free in the self-gravitating stellar case. The case of general baroclinic equilibrium states will be 
studied in Paper II. 

Results. Those theoretical results are applied to realistic stellar cases, namely to the solar radiative core and to the envelope of an Ap 
star, and discussed. In both cases we assume that the field is initially confined in the stellar radiation zone. 

Key words. Magnetohydrodynamics (MHD) - Plasmas - Magnetic fields - Sun: magnetic fields - Stars: magnetic fields 
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I— ^ 1. Introduction 



Spectropolarimetry is nowadays exploring the stellar magnetism 
acros s the whole H ertzspru ng-Russel diagram ([Donati et al. 



1997, 2006; Neiner 2007] |Landstreet et al.||2008t [Petit et al. 



2008). Furthermore, helioseismology and asteroseismology are 
providing new constra ints on internal transport processes occur- 
ing in stellar interiors ( Turck-Chieze & Talon|2008| Aerts et al. 
|2008| ). In this context, even if standard stellar models explain 
the main features of stellar evolution, it is now crucial to go be- 
yond this modelling to introduce dynamical processes such as 
magnetic field and rotation to investigate their eflTects on stellar 
structure and secular evolution ( Maeder & Meynet|2000 , Talon] 
[2008]). To achieve this aim, secular MHD transport equations 
have been derived in order to be introduced in stellar evolu- 
tion codes. They take into account in a coherent way the inter- 
action between diflTerential ro tation, turbulence, meridional cir- 
culation, and magnetic field ( |Spruit|[2002[ Maeder & Meynet 



Lydon & Sofia] 1 995 [[Couvidat et al.|20Q3||Li et al.|2006t[Duez 
et al.|20081|Li et al.|2009t|Duez et al.|2010| ). 

However, an infinity of possible magnetic configurations 
can be investigated since the diflTerent observation tech- 
niques only lead to indirect indications on the internal field 
topologies through the surface field properties they provide. 
Furthermore, since the simplest geometrical configurations like 
purely poloidal and purely toroidal fields are known to be unsta- 
ble (|Acheson|T978l|Tayler|1973, Markey & Tayler, 1 97^ [T9741 



Goossens & Veugelen|1978[|Goossens & Tayler] 1980[|Goossens 



et al. 



2ooa 



1981 



van Assche et al.||1982[ |Spruit||1999t [Braithwaite 

2007| ), the best candidates for stable geometries are mixed 

poloidal-toroidal fields ( Wrigh? p^73| [Markey & T^yler|p^74l 
|Tayler|1980t [Braithwaite 20 09 ). 



Therefore, it is necessary to go on the track of possible 
stable magnetic configurations in stellar interiors to evaluate 
their eflTects on stellar structure and to use them as potential 
initial conditions to study secular internal transport processes. 



^2004; Mathis & Z ahn||2005|), while non-linear numerical simu 
lations provide new insight on these mechanisms (Charbonneau 
& MacGregor[T993„Rudiger & Kitchatinov|1997„Garaud.2002, 



Brun & Zahn|2006| ). If we want to go further, the simplest mod- 
ifications of static structural properties such as density, gravity, 
pressure, temperature, and luminosity induced by the magnetic 
field have also to be systematic ally quantified as a function of the 
field geometry and strength ( |Moss|1973t[Mestel & Moss|1977| 



In this work, we thus revisit the pioneer works b y [Ferraro 
(1954); Mestel (1956); Prendergast (1956) and Woltjer[ ( [1960| 7 
^Ferraro, ( ,1954^ studied the equilibrium configurations of an 
incompressible sta r with a purely poloida l field. [Prendergast 
(1956) (see also Chandrasekhar T956a| [Chan drasekha r & 
Prendergast 1956; Chandrasekhar 1956b| then extended the 
model to take into account the toroidal field, by solving the 
magneto-hydrostatic equilibrium of incompressible spheres. The 



2 



V. Duez & S. Mathis: Relaxed Equilibrium Configurations to Model Fossil Fields 



obtained configurations seem to be relevant in regard of the most 
recent numerical simulations that may explain fossil fields in 
early- type st ars, white dwarfs or neutron stars (Braithwaite &| 
[Spruit 2004; [Braithwaite & Nordlu nd 2006; Braithwaite 2006) 
and of theoretical studies of their helicity relaxation ( Broderick| 
|& Narayan|2"Q08l|Mastrano & Melatosp OOS ). The main general- 
ization of Prendergast's work we achieve here consists in relax- 
ing the incompressible hypothesis in order to take into account 
the star's structure (see also |Woltjer||l960| ), which diff'ers as a 
function of its stellar type and of its evolution stage, and to de- 
rive the minimum energy equilibrium configuration for a given 
mass and helicity, which are then applied to realistic models of 
stellar interiors. 

Assuming that the Lorentz volumetric force is a perturba- 
tion compared with the gravity, we derive the non force-free 
magnetohydro static equilibrium. In this first article, we then 
focus on the barotropic equilibrium states family for which 
the possible field configurations and the stellar structure are ex- 
plicitely coupled; these may correspond to the numerical exper- 
iments by Braithwaite and collaborators. In this case, the prob- 
lem reduces to a Grad-Shafranov-like equation ([Grad & Rubin] 
[T958l|Shafranov|1966t[Kutvitskii & Solov'ev|1994| ), similar to 
the one intensively used in fusion plasma physics. We then fo- 
cus on its minimum energy eigenmodes for a given mass and 
helicity, which are derived and applied in order to model re- 
laxed stellar fossil magnetic fields which are found to be non 
force-free. Arguments in favor of t he stability of the obtained 
configurations are finally discussed ( Wright fl973[ Tayler'198()| 
jBraithwaite 2009 ; Reisenegger 2009) and we compare their 
properties with those of relaxed fields obtained in numerical sim- 
ulations ( Braithwaite 2008| ). The case of ge neral baroclinic equi- 
librium states will be studied in Paper II (Wright 1969[ Moss 
[1975). 



2. 1 . Magnetic field configuration and the 
magnetohydrostatic equilibrium 

If we consider only the axisymmetric case, where all physical 
variables are independent of the azimuthal angle (^), B (r, 6) can 
be written in the form 



B = (r, ^) X + F(r,0)e^, 

rsmO 



r sm 



(2) 



which is divergenceless. ^ and F are respectively the poloidal 
flux function and the toroidal potential; (r,6,(fi) are the usual 
spherical coordinates and {^k}k=r,e,(f their unit- vector basis. 
Finally, the poloidal component of the magnetic field (Bp) is 
such that Bp • = 0, so it belongs to iso-^ surfaces. The mag- 
netohydrostatic equilibrium (Eq.[T]) implies that the poloidal part 
(in the meridional plane in the axisymmetric case) of the Lorentz 
force (/^Xp) balances the pressure gradient and the gravitational 
force, which are also purely poloidal vectors, while, in the ab- 
sence of any other force, its toroidal component (/^£^ = i7:^e<^) 
vanishes. We thus have 



Using Eq. ([2]), we obtain 
1 



(3) 



fiQr^ sin^ 6 



1 



+ -(F5,F + 5,^A*^)e, , 



with dx = d/dx, and 



sin 6 



(4) 



(5) 



2, The Non Force-Free Magneto-Hydrostatic 
Equilibrium 

In this work, we focus on the magnetic equilibrium of a self- 
gravitating spherical shell to model fossil fields in stellar interi- 
ors. To achieve this goal, we start from 



= -VP-pVV + Fx, where Fj:=jx B. 



Eq. ([T]) must be satisfied in the interior of an infinitely conduct- 
ing mass of fluid in the presence of a large-scale field together 
with the Poisson equation, V^V = 4 n G p, and the Maxwell 
equations, V -B = (Maxwell flux) and VxB = /doj (Maxwell- 
Ampere). 

P, p and V are respectively the pressure, the density and the grav- 
itational potential of the considered plasma. B is the magnetic 
field and j is the associated current, which is given in the classi- 
cal MHD approximation by the Maxwell-Ampere's equation, jiq 
is the magnetic permeability of the plasma and F£ is the Lorentz 
force. 



^ Barotropic states are such that their density and pressure gradi- 
ents are aligned. They can be convectively stable or not, depending on 
their entropy stratification. We will precisely introduce their definition 
in §2.2. 



is the usual Grad-Shafranov operator in spherical coordinates. 
On the other hand, since F£^ = 0, we get dr^dgF - do^dyF = 0; 
the non-trivial values for F are therefore obtained by setting 



F(r, 6) = F(^). 



(6) 



Then, we obtain FdrF = Fd^iFdr^ and FdeF = Fdx^Fde^, 
leading to the final expansion of the Lorentz force 



(1) Fx = ^(r,^)V^, 



where 



1 



yUor^ sin 6 



(Fd^F -h A*^) . 



(7) 



(8) 



Therefore, the poloidal component of Fr is nonzero a priori, the 
field being thus non force-free in this case. 

This point has here to be discussed. First 



Reisenegger 



( [2009 ) demonstrated that the magnetic field can not be force- 
free everywhere in stellar interiors (see the demonstration in 
the appendix A of his paper). Note i n this context that the 
"force-free" configurations obtained by [Broderick & Narayan 
(j2008j) verify this theorem because they have current sheets 
with a non-zero Lorentz force on the stellar surface. Moreover, 
Shulyak et al. ( 2007 2009") showed how the atmosphere of a 
CP star can be the host of a non-zero Lorentz force. Therefore, 
from now on, we consider the non force-free equilibrium. 
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Now, if we take the curl of Eq. ([T]), we get the static vorticity 
equation 



VpxVP _ iFr 
-J—- = V X ' 



(9) 



that governs t he ba lance between the barocHnic torque (left- 
hand side; see |Rieuto rd ( 2006 ) for a detailed description) and 
the ma g netic source term. Then, as has been emphasized by 
Mestel (1956]), the different structural quantities such as the 
density, the gravitational potential, and the pressure relax in 
order to verify Eq. ^ for a given field configura tion (see [Sweet 
( [T950l ); possl ( |1975| ), and [Mathis & Zahnl ( |2nQ5l ) §5.). Thus, the 
choice for ^ is left free. 



2.2. The barotropic equHibrium state family 

Magnetic initial configurations are one of the crucial unan- 
swered question for the modelling of MHD transport processes 
in stellar interiors. To examine this question, Braithwaite 
and collaborato rs ( [Braithwaite & Spruit||2004( [Braithwaite & 
Nordlund|2006| ) studied the relaxation of an initially stochastic 
field in models of convectively stable stellar radiation zones. 
The field is found to relax, after several Alfven times, to a mixed 
poloidal-toroidal equilibrium configuration, which then diflTuses 
towards the exterior. 

We choose here, using an analytical approach, to find such 
field geometries, which are governed at the beginning by the 
magnetohydrostatic equilibrium. 

To achieve this aim, we focus in this first article on the partic- 
ular barotropic equilibrium states (in the hydrodynamic meaning 
of the term) for which the field configuration is explicitely cou- 
pled with the stellar structure, since we have in this case 



Vpx VP 



Vx 



0. 



(10) 



Those are the generalization of the Prendergast's equilibria 
which take into account the compressibility and which have been 



studied in poly tropic cases by Wo ltj er j (l 960| ) ; [^ntzel[ ( 1 96 1 1 ) ; 
Roxburg h (1966); |Monagh an (1976}r~^ 

Let us first recall the definition of the barotropic states. In 
fluid mechanics, a fluid is said to be in a barotropic s tate if the 
following condition is satisfied (see Pedlosky ( 1998| in a geo- 
physical context and ,Zahn,(1992j in a stellar one): 

VpxVP = 0; (11) 

in other words, the barocHnic torque in the vorticity equation 
(Eq. [9]) vanishes. Then, the surfaces of equal density coincide 
with the isobars since the density and the pressure gradients 
are aligned. This does not imply any question of equation of 
state, which in stellar interiors can take the most general form 
P = f (p,T,"') (T being the temperature). Moreover, this does 
not presume anything about the stratification of the fluid, which 
can be stably stratified or not. For example, a star in solid body 
rotation is i n a barotrop ic state (as opposed to barocHnic) (see 
once again |Zahn 1992). Let us illustrate this point with the 
simplest case of a non-rotating and non-magnetic stellar radi- 
ation zone. In this case, the hydrostatic balance is given by 
VP/p = -Vy = g. If we take the curl of this equation, we obtain 
the stationary version of the thermal- wind equation: 



and the star is thus in a barotropic state in the hydrodynamic 
meaning of the term. 

This has to be distinguished from the point of view of ther- 
modynamics where a barotropic equation of state is such that 
P = f(p) while a non-barotropic equation of state is such that 
P = f(p,Tr"). 

Then, it is clear that a fluid with a barotropic equation of 
state is automatically in an hydrodynamical barotropic state. 
However, in the case of a fluid with a non-barotropic equation 
of state, the situation is more subtle: in the case where the curl 
of the volumetric perturbing force vanishes (i.e. V x (F^/p) = 0) 
the fluid is in an hydrodynamical barotropic state, while in the 
general case it is in a barocHnic situation. Then, a fluid with 
a non-barotropic equation of state can be in a barotropic state 
even if it is only for a specific form of the perturbing force. 
In this first work, we choose to examine the first equilibrium 
family in which the Lorentz force verify the barotropic balance 
described by Eq. ( pT) in a stably stratified radiation zone. The 
second general case (cf. Mestel 1956) will be studied in paper II. 

Stellar interiors, except just under the surface, are in a regime 
where jS = P/Puag ^ 1, ^Mag = B^K^J^o) being the plasma's 
magnetic pressure. On the other hand, in the domain of fields 
amplitudes relevant for classical stars (i.e. the non-compact ob- 
jects), the ratio of the volumetric Lorentz force by the gravity 
is very weak. Therefore, the stellar structure modifications in- 
duced by the field can be considered as perturbations o nly from 
a spherically symmetric background ( [Haskell et al.|2008j ). Then, 
we can write p ^p-\-p, where p and p are respectively the mean 
density on an isobar, which is given at the first order by the stan- 
dard non-magnetic radial density profile of the considered star, 
and its magnetic-induced perturbation on the isobar (withp ^p). 
Thus to the first order, Eq. ([TO]) on an isobar becomes 



P \ P 



(13) 



where the eflfective gravity (g^ff), such that VP = pg^^, has been 
introduced. This gives, using Eq. ^ 



which projects only along e<^ as 

SO that there exists a function G of ^ such that 
P 

Then, Eq. ^ leads to the following one ruling Y 
A*'¥ + F 5vp [F CP)] = -juor^ sin^ £(p G . 



(14) 



(15) 



(16) 



(17) 



This equation is similar to the well-known Grad-Shafranov 
equatioij^ which is used to find equilibria in magnetically con- 
fined plasmas suc h as those in tokamaks or in spheromaks ( GradJ 
& Rubin 1958; Sha franov[T 966). However, here the source term 
is diflTerent and is directly related to the internal structure of the 



VpxVP 



-V X [VV] = 



(12) 



^ The usual Grad-Shafranov equation is given by: 
A*^ + F C¥) [F (^)] = -hqv^ sin^ 68^ [P (^)], 
where the pressure P is prescribed in function of ^. These describes the 
equilibrium between the magnetic force and the pressure gradient only. 
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star through its density profile (p) (the general form of the Grad- 
Shafranov equation in an astrophysical context is discussed for 
example in Heinemann & Olbert (1978 ) and Q gilvie| ( [1997| )). 



Moreover, since the field has to be non force-free in stellar inte- 
riors G ^ (see the previous discussion in §2.1. and Eqs.[T6j|7] 
andjS]). 

It is only applicable to the case of the barotropic state family. 
The equations for the general case will be studied in paper II. 

2.3. F and G expansion 

Let us now focus on the respective expansion of F and G as a 
function of ^. 

First, since F is a regular function, we can expand it in power 
series in ^: 



(18) 



the Xi being the expansion coefl&cients that have to be deter- 
mined and R a characteristic radius which will be identified 
below. On the other hand, must be regular at the center of 
the sphere; the first term (/ = 0) of the previous expansion is 
then excluded (cf. Eq. [2]), the above expansion thus reducing to 

In the same way, G can be expanded as 



G(^) = YjPj^K 



Then, Eq. ( 17 ) becomes: 



A*^ + Yj = -^or' sin^ Op £ A^^, 



(19) 



(20) 



k>0 



where A,^ = Z/i>o Z/2>o {^'2^/i^/2^/i+/2-i,^}' ^ being the usual 
Kronecker symbol. This is the generalization of the Grad- 
Shafranov-type equation obtained by Prendergast (1956) for the 
barotropic compressible states. 

Thus, having assumed the non force-free barotropic 
magneto-hydrostatic equilibrium state leads to undetermined ar- 
bitrary functions (F and G) that must be constrained. To achieve 
thi s aim, we follow the method given in the axisymmetriccase 
by [Chandrasekhar & Prendergast] ([1958 ) and Woltjer] ( |1959b| ) 
that allows to find the equilibrium state of lowest energy com- 
patible with the constancy of given invariants for the studied ax- 
i symmetric system. 



3. Self-Gravitating Relaxation States 

3.1. Definitions and axisymmetric invariants 

We first introduce the cylindrical coordinates (s, (p, z) where s = 
r sin 6 and z = r cos 6. Then, B given in Eq. (|2]) becomes 



B{s,z)= -V^(^,z)x ^^ + -F{s,z) e^. 
s s 



(21) 



Then, we define the potential vector A {s, z) = (s, z) e<^ such 
that Bp = V xA and we get 



F ^ 
B = V xA -h — where (s,z) = — . 

s s 



(22) 



Next, we insert the expansion for the ma gnetic field B use d by 
Chandrasekhar & Prendergast] ( p^58l ) and |Woltjer| ( |1959a|b| ): 



B = -s d,^ (s, z) Cs + [s^^ (s, z)] h + sT (^, z) (23) 
where ^k}k^s,ip,z cylindrical unit- vector basis and where we 

(24) 



identify using Eq. (21 
^ = s^^ and F = s^T. 



The Grad-Shafranov operator applied to ^ can then be expressed 
as follow 











dss--ds+dzz 






s 





where A5 = dss + -ds + ^zz- 



(25) 



We now introduce the two general families of invariants 
of the barotropic axisymmetric magneto-hydrostatic equilibrium 
states, which have been int roduced by Woltjer, ( (1959bj for the 
compressible case (see also Wentzel|196Q| ): 



Il,n = j^^n (^'O) pd^= {s^^f p dT, 



(26) 



(27) 



where M„ and A^^ are arbitrary functions that have to be speci- 
fied. Those latters are conserved as long as 



B-7r = (i.e. (D = r = 0) 
on the boundaries. 



(28) 




Fig. 1. Schematic representation of the two coordinates systems 
used and of a constant ^ surface. The invariants of the axisym- 
metric system are the total mass of the considered stellar radia- 
tive region (Mrz), the mass enclosed in a constant ^ surface, 
the toroidal flux (^7^) associated with the toroidal magnetic field 
(B^), and the global hehcity (*K). 
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3.2. Fossil fields barotropic relaxation states 

Let us first concentrate on Ij.q and A^^ relevant for fossil fields 
relaxation. First, if we set A^o (^s^^^ = 1, we obtain 



= r Td'V = 2n f B^dsdz 

JS 



(29) 



that corresponds to the conservation of the flux of the azimuthal 
field across the meridional plane of the star {T^p) in perfect 
axisymmetric MHD equilibria. 



Then, if we set Ni i^s^ = O, we get 



(30) 



where we thus identify the magnetic helicity (TY; see §5.1.) of 
the field configuration which is a global quantity integrated over 
the volume of the studied radiation zone. 

Let us briefly discuss the peculiar role of this quantity in 
the search of stable equilibria. As emphasized by Spruit ( 2008 ), 
the magnetic helicity is a conserved quantity in a perfectly con- 
ducting fluid with fixed boundary conditions. However, in re- 
alistic conditions, rapid reconnection can take place even at 
very high conductivity, especially when the field is dynami- 
cally evolving, for example during its initial relaxation phase. 
Nevertheless, in laboratory experiments, like for example in 
spheromaks, the helicity is often observed to be approximately 
conserved, which leads to the existence of stable equilibrium 
configurations. In fact, if the helicity is conserved a dynamical 
or unstable field with a finite initial helicity (*Ho) cannot decay 
completely, the helicity of vanishing field being zero. This is 
precisely what has been observed in the nu merical experiment 
performed by [Braithwaite & Spruit] ( |2004| ) and Braithwai te &| 
[Nordlund"^ ^2006), where an initial stochastic field with a finite 
helicity decays initially but relaxes into a stable equilibrium. 

In the context of laboratory l ow-jS plasmas, this process has 
been identified by Taylor (1974]) and is thus called the Taylor's 
relaxation. 



For this reason, we now search as [Chandrasekhar &| 
Prendergast (1958) the final state of equilibrium, which is the 
state of lowest energy that the compressible star, preserving its 
axisymmetry, can attain while conserving the invariants Ji;„, 
■^1,0 = ^ and Ij-i = 'H/2 in barotropic states. 

To achieve this, we thus introduce the total energy of the 
system 



-^p[V-^2^(]\d^ 



" \J [-^^OAsO s^T^] + p [y + 2^(]^ d^. 



(31) 



where U is the specific internal energy p er unit mass ( [Woltjer 
T9581 |1959bt [Broderick & Narayan|[jD58] ). To obtain the mini- 
mal energy equilibrium state for the given invariants Ij-n and a 
given helicity and azimuthal flux, we thus minimize E with re- 
spect to Ii-n, J^i^o and I^-i . This leads, introducing the associated 



Lagrangian multipliers {ai-n,aj[-o,aj[-i), to the following condi- 
tion for a stationary energy: 



+ ^ ai-n Sli-n + ^ Cll;q ^Il;q = 0. 
n q^O 



(32) 



Following the method describ e d in jChandrasekhar & 

Prendergast] ( [19581 ) and jWoltjer] ( |1959b| ), we express 6E 
and 6Ij-r in function of ^O, 6T and 6p. Since these variations 
are independent and arbitrary, their coeflScients in the integrand 
of Eq. ([32]) must separately vanish, which gives 



Po 



dMn{s^^) 1 dN,{s^^) 



_ dMn{s^^) 
Ph^'-^^ d(.^cD) 



^ 1 



Po 



(33) 



(34) 



^=0 



These equations thus describe the minimal non force-free energy 
equilibrium states for a given helicity and azimuthal flux. From 
Eq. ( [Tt] ), we identify that F Q¥) is now constrained while G Q¥) 
required to ensure the non force-free character of the field is still 
arbitrary. 

Let us now consider the first invariants family (Ji) given in 
Eq. ( [26] ) which are thus necessary to constrain G First, the 
non-magnetic global quantity, which is an invariant of the con- 
sidered equilibrium, is the total mass of the stellar radiation zone 
Mrz. We thus set Mq (^^O^ = 1, leading naturally to consider the 
mass 



I 



i;0 



r pdT 



Ml 



RZ- 



(35) 



However, since dMo (^s-^O^ d (^-^O^ = 0, we have thus to con- 
sider the highest-order invariant 



I 



pdT 



(36) 



because of the non force-free be haviour of the field . This last 
invariant has been introduced by Prendergast ( 1956*) in the ax- 
isymmetric non force-free magneto-hydrostatic incompressible 
equilibrium and it corresponds to the mass conservation in each 
flux tube described by the closed magnetic surface s^(^ = c^^^. 

Furthermore, the considered radiation zone is stably strat- 
ified. Since in stellar interiors the magnetic pressure is very 
much less than the thermal one, the Lorentz force have only a 
negligible eflTect on the gas pressure (J3^l). Moreover, energy 
is required to move fluid elements in the radial direction because 
work has to be done against the buoyant restoring force which 
is thus very strong compared to the magnetic one. Therefore, 
the radial component of the displacement (^) which takes place 
during the adjustment to equilibrium is inhibited ^ • 'e^ ~ and 
V • (pf ) ^ due to the anelastic approximation justified in stellar 
radiation regions. Therefore, as emphasized by jBraithwaite 
(2008 ), the mass transport in the radial direction is frozen (no 
matter can leave or enter in the flux tube) and Ij-i can be used 
as a supplementary constrain in our variational method. 
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From Eqs. ( 33p4 ), we thus get the following equations de- 
scribing the barotropic axisymmetric equilibrium state of lowest 
energy that the compressible star can reach while conserving its 
radiation zone mass (Ji;o = Mrz), the mass in each flux tube 
(Ji i), the flux of the toroidal field (Ji o = T^p), and a given he- 
licity (Jn;i = ^/2): 



— A5O = ai-i p + a^-i 
— r = -a]i;iO -. 



(37) 



(38) 



Since the azimuthal field is regular at the origin, we get from Eq. 
(23 ) aj[-o = 0. Eliminating T between Eqs. ( 37p8 ), we obtain 



(39) 



that becomes, multiplying it by and using Eq. (24 ) & (25 ) 
A*^ + [Moaj-if ^ = yuo ai-i p ? sin^ 6. (40) 



We thus identify in Eq. ( 20 ) 



j=0 \a,,=-j-^AJR^ 



(41) 



where we have constrained the initial arbitrary functions of the 
magnetohydrostatic equilibrium 



F (^) = -jioajL-i^ and G (^) = -ai-i . 
So, it reduces to 



(42) 



(43) 



the values of the real coefl&cients Ai and j3o being thus controled 
by the helicity (*K) and the mass conservation in each axisym- 
metric flux tube defined by ^ = c^^^ because of the non force-free 
stably stratified behaviour of the reached equilibrium. 

This corresponds, as has already been emphasi zed, to the 
lowest energy equilibrium state for a given helicity ( B ellan|2000[ 
Broderick & Narayan|2008 ). The equilibrium state ruled by Eq. 



(|43|) are thus the generalization of the Taylor relaxation states 
in a self-gravitating star where the field is not force-free (i.e. 
V X B ^ a B). Note also that some non force-free relaxed states 



have been id entified in plasma physics (Montgomery & Phillips 
p?88, 1989 i Dasgupta et al.|2002t [Shaikh et al.|2008| ) that should 
be studied in a stellar context in a near future. 

Let us note that in the case where is not con sidered 
(ai l = y^o = 0) we recover thej Chandrasekhar ' ('1956a') force 
~ ?2 



Using Green's function method ( [Morse & Feshbach|1953 ) , we 
then obtain the particular solution associated with S: 



% (r,x) = -yUoySoX^ 



2/-h3 



^ .Jsup [2(/+l)(/ + 2) 



sup 



x[l-x^)c^/^ (jc), 
where 

J^/(^) = j|<S(^,/)q'/'(x)d/; 



R, 



sup 



(46) 



(47) 



ji and yi are respectively the spherical Bessel functions of the 
first and the second kind (also called Neumann fu nctions) while 
C^^are the Gegenbauer polynomials ( ^Abramowitz & Stegun 
'1972[). Rm and ^sup, which are respectively the bottom and the 



top radius of the considered radiation zone, are introduced and 
we identify R = R^up- 

These functions (7/, yi, and C^/^) are respectively the radial 
and the latitudinal eigenfunctions of the homogeneous equation 
associated with Eq. ( [44| ) : 



0. 



(48) 



Then, if we express the solutions of this equation as 
Z/ // (r) 81 (0), we get respectively 



and 

dr2 



dx2 



+ (/+l)(/ + 2)g, = 



sup 



(/+!)(/ + 2) 



fi = 0. 



giving 

gi = {l-^)c^'Hx) 
and 



(49) 



(50) 



(51) 



sup 



sup 



^sup \ ^sup / 



(52) 



K[ and being real constants. A[ are eigenvalues that allow to 



free limit (see also Marsh| |1992[ for a generalization of the ^^^^^^ the' boundary conditions discussed hereafter. One has to 



solutions) and the usual Taylor's states for low-/? plasmas. The 
Prendergast's model is recovered assuming a constant density 
profile (incompressible). 



3.3. Green's function solution 

We are now ready to solve Eq. ( [43] ). If we introduce x = cos 6 
and if we set S (r, 6) = -jiq JSq p P~sin^ 6, it is recast in 



notice that has to vanish in order to preserve the regularity of 



the solution at the center. Applying this to Eq. (43), we finally 
obtain the general solution 



^(r,6>) 



sm 



/=0 



where 



drr + 



1 
7?2 



(44) 



(45) 



^sup \ ^sup / Jr L \ ^sup / 
rtsup \ «sup/jRi„fL \ ^sup/ 



(53) 
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Notice that the particular solution for the poloidal flux 
function (^p^ presents a dipolar geometry, owing to its an- 
gular dependence that follows the one from the source term 
S = -yuo y^o P sin^ 6. 

Neglecting the density, we end up with the linear homogeneous 
equation, whose solutions are Cha ndrasekhar-Kendall functions 
( [Chandrasekhar & Kendall] |1957| ). Moreover, the source term 
being constituted only of a dipolar component, all the non- 
dipolar contributions are, according to this model, force-free. 



The magnetic field is then given for r < 7?sup by 



B 



1 



sin 6 



rsinO 



Rsup r sin ^ 



After a few manipulations, we can then express the current den- 
sity as 



Jp 



-VX^T 



Bp 



/V = — V x^p 




= — -^T+y^oP^sin^e,^, 



(55) 



(56) 



non force-free 



^ (7? = ,0) = 0, which gives the two independent equations 
for / = 



3 A 



Rn 



and 



(58) 



(59) 



we here focus on the dipolar mode which is known to be the 
lowest energy per helicity ratio state (cf. Bro derick & Narayan 
12008 ). These can be formulated so that one first determines the 



(54) value of /Ij ' according to 



71 



Rc 

iO,( ^ 



d^ = 



(60) 



and next computes following (59 ). 



In the case where there is no convective core, as for example 
in central radiation zones of late-type stars such as the Sun, Eqs. 
(59)) & ([60]) must be appHed setting R^ = 0. 



where we recognize in the first term of the right hand side the 
force-free contributions and in the second the non force-free one, 
fully contained in the toroidal component. 

The Lorentz force can as a matter of fact be written in the 
very simple form: 



3.4. Configurations 



(57) 



The boundary conditions for ^ which determine possible values 
for K[ and /l^'^ have now to be discussed. Two major types of 
geometry are relevant for large-scale fossil magnetic fields in 
stellar interiors: initially confined and open configurations. 



3.4.2. Open configurations 

This corresponds to the case of fields that match at the stellar 
surface (at r = R^, R^ being the star's radius) with a potential 
field as observed now in some cases of early-type stars such as 
Ap stars. Then, we have ^ext = VOm, being the associated 
potential. 

In the case studied here, we focus on the first configuration 
(initially confined) since the search of relaxed solutions for given 
i^i;0, ^]i;0, and I J- 1 assumes that 5 • Cr = on the stellar 
radiation zones boundaries. This initial confined configuration 
will then become open one through Ohmic diffusion as in the 
Braithwaite and collaborators scenario. 



3.4.1 . Initially confined configurations 

Let us first concentrate on the simplest mathematical solution 
in the case of a central radiation zone that initially cancels ^ 
both at the center (/?inf = 0) and at a given confinement radius 
(^sup = Rc)- 

Then, if we choose to cancel the K[ coefl&cients 
for every /, the condition ^(0,6) = is verified since 

lim^^o ^ ji ^^i'^ ^ ) - ^- However, if we look at the magnetic 

field radial component behaviour at the center, it is easily shown, 

using Eq. ( 54 ), that if = it is given by lim^^o r~ ^ ji (a^/ ^ j 

which does not cancel so that 3^(0,6) ^ = C cos 6 (where 
C G R*). Therefore, this solution is multivaluated, thus physi- 
cally inadmissible and + 0. 

Then, we consider the general case of a field initially con- 
fined between two radii Ri^f - Rc^ and R^^p = Rc^ , owing to the 
presence of both a convective core and a convective envelope (as 
it is the case e.g. in A-type stars). We impose ^ (7?^ , ^) = and 



4. Application to Realistic Stellar Interiors 

To illustrate our purpose, we apply our analytical results (i) to 
model an initial fossil field buried below the convective enve- 
lope of the young Sun on the ZAMS, and then (ii) to model an 
initial field present in the radiation zone of a ZAMS 2.40 M© 
magnetic Ap-star, whose lower and upper radiation-convection 
interfaces are respectively located at 7?^ = 0.111 R^ and at 
Rc^ = 0.992 R^. In the first case, the parameter fio is deter- 
mined such that the maximum field strength reaches the ampli- 
tude of ^0 = 2.1 MG, which is one of the upper limits given 
by F riedland & Gruzinov| ( |2004) for the present Sun's radia- 
tive core. In the second case, it is obtained such that it reaches 
the arbitrary value of = lOkG. This value is approxima- 
tively of the same order of magnitude that the mean surface am- 
plitude observed using spectropolarimetry for magnetic Ap-star 
which exhibits strong external dipolar magnetic behaviour (such 
as HD12288 ( Wade et al.|2QQ0| )). We thus assume implicitly that 
such an initial confined internal field is a potential prelude to the 
multipolar one observed now at the surface, the latter state be- 
ing acheived after a diflfusive process that will be studied in a 
forthcoming paper. 
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4.1. Fossil fields buried in late-type stars radiative cores 

The young Sun model used as a reference is a Cesam non- 
rotating standard one (Morel 1997 |, following inputs from the 
work of |Couvidat et al.,(2QQ3) and |Turck-Chieze et al.| ( |2QQ4| ). 

In Fig. |2] three possible configurations for ^ are given. We 
choose those corresponding to the first, the third and the fifth 
eigenvalues (given in table [T]). Those are the generalization of 
the well-known Grad-Shafranov equation linear eigenmodes ob- 
tained in the force-free case (cf. Marsh 1992). The field is then 
of mixed-type (B^p is given for /l^'^), both poloidal and toroidal, 
and non force-free, properties already obtained by Prendergast 
( (1956j ) in the incompressible case. The respective amplitudes ra- 
tio between the poloidal and the toroidal components will be de- 
scribed in §5., where the possible stability of such configurations 
will be discussed. 



Table 1. Eigenvalues of the first five equilibria for the two con- 
figurations illustrated. 



Eigenvalue 


Solar case 


Ap star case 




5.276 


4.826 


f 

4' 


9.157 


8.657 


12.951 


12.444 


16.290 


16.174 


19.839 


19.849 



4.2. Fossil fields in early-type stars 

Respective corresponding possible configurations in the case of 
an Ap star are given in Fig. [3] 

The model is typical of an A2p-type star, with an initial 
mass Ma = 2.40 M©. The solar metallicity is chosen as the 
initial one and the model is taken on the ZAMS, its luminosity 
being = 38.0 L©. 

Obtained configurations are then mixed poloidal-toroidal 
(tw isted) fields which may be stable in stellar radiation zone s 
(cf . [Braithwaite & Spruit|2004l [Braithwaite & Nordlund|2QQ6i ). 
Their configurations are thus given in both cases by concentric 
torus, the neutral poloidal points (where (9^^ = dg^ = so that 
Br = Bq = 0) positions being function of the internal density 
profile of the star. 

Let us emphasize here that the original approach of this 
work first consists in deriving the Grad-Shafranov-like equation 
adapted to treat the barotropic magnetohydrostatic equilibrium 
states for realistic models of stellar interiors. Such approach 
has already been applied to investigate the internal magnetic 
configurations in polytropes and in compact objects such that 



& Melatos 20041 Tomimura & Eriguchi||2005[ |Yoshida et al. 


2006[ Haskell et al.||2008 


Akgun & Wasserman||2008 


Kiuchi 



Then, the obtained arbitrary functions are constrained with 
deriving minimal energy equilibrium configurations for a given 
conserved mass, azimuthal flux and helicity that generalizes the 
relaxation Taylor's states to the self-gravitating case where the 
field is non-force free. 



5. Links between the field's helicity, topology and 
energy 

5.1. Helicity vs. mixity 

Let us express the magnetic field (B) in terms of magnetic stream 
functions (for the poloidal component of the field) and (for 
its toroidal part): 



5 = V X [V X [^P (r, e)er] + (r, ^) e,] . 



(61) 



Next, the vector potential A is given by the relation B = V x A. 
Knowing that in the confined case the gauge choice is inconse- 
quential, we can identify without further ado 



A = Vx[^p(r,^)'5;]+^T(r,^)e,. 



(62) 



The magnetic stream functions are then projected on the spheri- 
cal harmonics 



^Ar,e) = ^^^(r)y,^^). 



(63) 
(64) 



From now on, we use Einstein summation convention where 
B{ = YaI^^ ^nd the vectorial spherical harmonics basis 
(r^ {6) , {6) , {6)) such that any axisymmetric vector field 
M(r, 6) can be expanded as 



M(r, 6) = <(r) (6) + v'^(r) (6) + w^(r) (6) , 



(65) 



where the vectorial spherical harmonics R^ (0), (6), and (0) 
are defined by: 

R? (6) = (6) (6) = VsY^ (0) , (0) = x R^ (6), (66) 



the horizontal gradient being defined as = eedg. (cf. Rieutord 
( 1987) ). Since 



V X (^p?,) = V X (^^ Y^er) = V X (^^ R?) 
we get from Eq. ( [62] ) 



^0 rpO 



(67) 



(68) 



On the other hand, we have 



^0 



r r 

from which we finally obtain the expression for the helicity 



(69) 



r (r + 1) 



\&a?dr. (70) 



At this point, we define a "poloidal helicity" defined by 



Ap . 5p d^V 



and a "toroidal helicity" by 

X 



AsB. d'V. 



(71) 



(72) 
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^, ■^'/^'o isocontours 




B, (G) 




■ 


2.00E + 06 1 _ 






l,B0E + 06 






1.60E + 06 






0.5 - 






1.40E + 06 






1.20E + 06 / 






' 






l.OOE + 06 ^ 






8.00E + 05 \ 






6.00E + 05 

0.5 - 




1 


4.00E + 05 






2.00E + 05 






O.OOE + 00 1 



7 isocontours 



■ 



l.OOE + 00 

8.00E-01 

6.00E-01 

4.00E-01 

2.00E-01 

O.OOE + 00 ^ 

-8.25E-02 

-1.65E-01 

-2.47E-01 

-3.30E-01 

-4.12E-01 




7 values 



I 



.OOE + 00 
8.00E-01 
6.00E-01 
4.00E-01 
2.00E-01 
O.OOE + 00 
-2.00E-01 
-4.00E-01 
-6.00E-01 
-8.00E-01 
-l.OOE + 00 



I 



I 



l.OOE + 00 
8.00E-01 
6.00E-01 
4.00E-01 
2.00E-01 
O.OOE + 00 
-4.00E-02 
-7.99E-02 
-1.20E-01 
-1.60E-01 
-2.00E-01 



Fig. 2. Upper panel: (left) toroidal magnetic field strength in colorscale and normalized isocontours of the poloidal flux function ^ in 
meridional cut in the solar case where the field is buried in the radiative core (below 0.1 26R^) for the first equilibrium configuration 



(/l^'^); (right) anisotropy factor y (Eq. 88|). Lower panel: normalized isocontours of the flux function ^ in meridional cut (left) for 



the third possible eigenvalue (A^{^) in the same case, (right) for the fifth possible eigenvalue (A^{^). The dashed circles indicate the 
radiation-convection limits. 



^, ■^'/^'o isocontours 




B. (G) 






l.OOE + 04 1 




1 


9.00E + 03 






8.00E + 03 

0.5 - 






7.00E + 03 






6.00E + 03 1 






5.00E + 03 \ ol 






4.00E + 03 ' 






3.00E + 03 




1 


0.5 - 

2.00E + 03 






l.OOE + 03 






O.OOE + 00 1 



7 isocontours 



7 values 









l.OOE + 00 ^ _ 




1 


8.00E-01 






6.00E-01 

0.5 - 

4.00E-01 






2.00E-01 (i 






O.OOE + 00 \, 






-1.73E-01 ^ 






-3.46E-01 

0.5 - 






-5.19E-01 






-6.93E-01 






-8.66E-01 1 




l.OOE + 00 
8.00E-01 
6.00E-01 

- 4.00E-01 

- 2.00E-01 

- O.OOE + 00 

- -1.43E-01 

I-2.85E-01 
-4.28E-01 
-5.71E-01 
-7.13E-01 



Fig. 3. Upper panel: (left) toroidal magnetic field strength in colorscale and normalized isocontours of the poloidal flux function ^ 
in meridional cut in the Ap star's case where the field is confined between 7?^ = 0. 1 1 1 and = 0.9927?* for the first equilibrium 
configuration (/l^'^); (right) anisotropy factor y (Eq. 88 ). Lower panel: normalized isocontours of the flux function ^ in meridional 



cut (left) for the third possible eigenvalue (A^{^) in the same case, (right) for the fifth possible eigenvalue (A^{^). The dashed circles 
indicate the radiation-convection limits. 



Since (r^, S^, T^) constitutes an orthogonal basis, 

r R^ S^dIl= r R^ T^dIl= r S^ X^d^l^O, 

JQ Jq. Jq. 



(73) -K 



we get from the previous expression: 
\£' {£' + l) 



Jo 



4 r n-K 

Jq. 
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and we verify that 
<H = -Hp + -Wt- 



(74) 



(75) 



From this expression (74) we can draw two conclusions: 



1 . A magnetic field has to be mixed (both poloidal and toroidal) 
to be helical; 

2. The poloidal and the toroidal helicities are equal. This can 
be verified, exploiting the orthogonality relations 



X 

and 

L 



(76) 



(77) 



(78) 



S£^£' being the usual Kronecker symbol. Then, we get 
Jo 

5.2. Helicityvs. energy 

Now, we focus back again on the helicity expression in terms 
of the poloidal flux function ^. The equations ( [55] ) and ([56]) are 
rewritten as 



R 



VX^T, 



R R 

= —yx[Bp-/uol3opr^ cos6>er]. 
We thus obtain the two vector potentials 
Ap = — (^p -fioPop cos e Cr) + VAp, 

At = -^^t + VAt, 



(79) 

(80) 

(81) 
(82) 



Ap and Aj being scalar gauge fields left free. When deriving 
the poloidal and toroidal helicities with the boundary condition 
^ • e"^ = at the surface, these ones disappear after integration 
and we find : 



(83) 
(84) 



So, introducing the poloidal and toroidal magnetic energies 

f/mag;P = jy ^ and t/mag;T = jy ^ (respectively), 
we obtain: 



IllqRI 1 



(85) 



where we identify Mxp = Jj-i = Jy p ^ d*V, and 



Qj -^mag;T. 



(86) 



Finally, adding these two last equations, we get the global re- 
lation between the helicity and the magnetic energy in the non 
force-free case 



(87) 



IuqRI 1 \ 

^=^r(^-ag-2/^0M^]' 

where we recognize in the second term the non force-free contri- 
bution, which is the first invariant: the mass enclosed in magnetic 
flux surface. 



5.3. Helicity vs. topology 

The / > 1 latitudinal modes contributions 

As shown by Broderick & Narayan (2008) for a set of modes 
/ ranging from 1 to 8 in the case of force-free solutions applied in 
an incompressible media, the first dipolar eigenvalue /l^'^ corre- 
sponds to the minimum energy configuration. Furthermore, from 
the Eq. (87 ), it arises directly that adding contributions from the 
higher multipolar components of the field (force-free) will result 
in adding a positive amount of magnetic energy to the total en- 
ergy, and this one will not be the minimal state. 



Lowest energy radial mode 

We plotted in Figure[4^ and[4]3 the poloidal, toroidal and to- 
tal helicity in the case of the Sun and of the Ap star. It clearly 
stems from this figure that the poloidal and the toroidal helicities 
are equal (cf . Eq. 78 ) for the eigenvalues given in Tab [T] (repre- 
sented by red diamonds). Moreover, the energy of the poloidal 
component of the field (t/mag;p) can be compared to the one cor- 
reponding to the toroidal part (t/mag;T) and we see that they are 
of the same order of magnitude. 

In Figure |5^ and|5|3 are represented the ratios fimag/'W for 
the poloidal, toroidal and global contributions, with and without 
the non force-free term, as a function of the parameter /l^'^ re- 
spectively in the case of the Sun and of the Ap type star. The first 
dipolar eigenvalue /l^'^ present the minimum energy compared 
with highest radial modes. It is thus the most probable config- 
uration achieved after relaxation and from now on we focus on it. 



6. Discussion 

6.1. Stability criteria 

First, it is interesting to examine the ratio of the field's poloidal 
component amplitude with its toroidal one. Then, we define the 
anisotropy factor (7) of the configuratioij^by 



r(^^) = 



Bp -h B^ 



where Bi 



+ Bi 



(88) 



It runs between -1 when the field is completely toroidal to 1 
when it is completely poloidal. In Figs.[2]&[3j it is shown for 

^ This can be inverted as 
^ = and ^ = IzZ, 

where the magnetic energy densities associated respectively with the 
poloidal field (Wmag;? = Bl/l/do) and with the toroidal one (Wmag;! = 
B\I2iiq) have been introduced. 
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Fig. 4. Normalized total, poloidal, and toroidal helicities as a 
function of the eigenvalue (/l^) in the case (a, top) of the young 
Sun and (b, bottom) of the studied Ap star. The red diamonds 
represent the eigenvalues C^'^) given in Tab 111 for which Eq. 
( [7^ is verified. Using Eq. ( [86] ), we directly deduce t/mag;T while 
1 7^g;P is given in purple. 

the first configurations obtained in the solar and in the Ap star 
cases. In both ones, the field is strongly toroidal (y ^ -1) in 
the center of the torus, which corresponds to the neutral point 
of the poloidal field (where we recall that = de^ = 0), 
while it is strongly poloidal (y ^ 1) around the magnetic axis 
of the star where the toroidal field is weak. Between those two 
regimes, both components have comparable strengths where 
y ^ 0. Then, proposed configurations may be stable since 
poloidal and to roidal fields can st abilize each other respectively 
(fright 1973) |Tayler| | 



stability analysis fof 
[Bernstein e'tari ( |1958 



T98Q| [Braith waite 2009). The complete 



owing the analytical method given in 
and using 3D numerical simulations will 



be achieved in a near future. 



6.2. Comparison to numerical simulations 

Next, let us compare in more details our analytical configuration 
to those obtaine d using numerical simulations (se e |Braithwaite 
& Spruit] ( [2004] ); [Braithwaite & Nordlund] ( [20061 ); iBraithwaite 
( I2008D ). 

Braithwaite and collaborators performed numerical magne- 
tohydrodynamical simulations of the relaxation of an initially 
random magnetic field in a stably stratified star. Then, this ini- 
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U„,g/Hp X fxji 
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Fig. 5. Magnetic energy/helicity ratios for the total, poloidal, and 
toroidal contributions as a function of the eigenvalue (/l^) with 
(first five curves) and without (last three curves) the non-force- 
free term in the case (a, top) of the young Sun and (b, bottom) of 
the studied Ap star. For the eigenvalues (A^/) represented by the 



red diamonds (cf. Tab[T]), Eqs. ( 78|87|85[86 ) are simultaneously 
verified. 



tial magnetic field is always found to relax on the Alfven time- 
scale into a stable magneto-hydrostatic equilibrium mixed con- 
figuration consisting of twisted flux tube(s). Two families are 
then identified: in the first one, the equilibria configurations are 
roughly axisymmetric with one flux tube forming a circle around 
the equator such as our configuration; in the second family, the 
relaxed fields are non-axisymmetric consisting of one or more 
flux tubes forming a complex structure with their axis lying at 
roughly constant depth under the surface of the star. Whether an 
axisymmetric or non-axisymmetric equilibrium forms depends 
on the initial condition chosen for the radial profile of the initial 
stochastic field strength \B\\ oc p^: a centrally concentrated one 
evolves into an axisymmetric equilibrium as in our configuration 
while a more spread-out field with a stronger conn ection to the 
atmosphere relaxes into a non-axisymmetric one. Braithwaite 
( 2008) indicates that, using an ideal-gas star modelled initially 
with a polytrope of index n = 3, the threshold is p ^ 1/2. 



Moreover, as shown in Fig. 7 in Braithwaite (2008 ), a selec- 
tive decay of the total helicity (*K) and of the magnetic energy 
(^Mag) occurs during the initial relaxation with a stronger de- 
crease of Uuag than that of "H. This hierarchy well known in 
plasma phy sics (see for example [Biskamp[ ( [1997| ) and [Shaikh 
et al. ( 2008| )) justifies the variational method used to derive our 
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configuration (Montgomery & Phillips ( 1988| )) while the intro- 
duction of 1 1- 1 is justified by the non force-free character of the 
field in stellar interiors ( Reisenegger 2009 ) and by the stratifica- 
tion which inhibits the transport of flux and mass in the radial 
direction (see §3.2. and Braithwaite|(j2008j)). 

Finally, note that our analytical configuration for which 
^^mag;p/f^mag ~ 0-45 verifies the stability criterion derived by 
[Braithwaite^ i^QQ9 ) for axisymmetric configurations: 



grav 



(89) 



mag 



where t/grav is the gravitational energy in the star and J{ a dimen- 
sionless factor whose value is of order 10 in a main-sequence star 
and of order 10^ in a neutron star while we expect 6^mag/^ grav < 
10"^ in a realistic star (see for example Duez et al. (2010)). 

Our analytical solution is thus similar to the axisymmetric 
non force-free relaxed solutions family obtain ed by Braithwaite^ 
& Spruit] ( [2QQ4| ) and |Braithwaite & Nordlund| ( |2006t . ' 



Those types of configurations can thus be relevant to model 
initial equilibrium conditions for evolutionary calculations 
involving large-scale fossil fields in stellar radiation zones. 
First, they can be used to initiate MHD rotational transport 
in dynamical stellar evolution codes where it is implemented 
(cf. Mathis & Zahn|[2005| ) . There, axisymmetric transport 
equations have been derived to study the secular dynamics of 
the mean axisymmetric component of the magnetic field, the 
magnetic ins tabilities being treated using phenomenolo gical 
prescriptions ( |Spruit|1999j [20021 [Maeder & Meynet|2 004 ) that 
have to be verified and improved by numerical experiments 
jNordlundl p006|) a nd subsequent works; |Zahn et al.| ([2007); 
jGellert et al.[ ( |2008) ). On the other hand, those can also be used 
as initial conditions for large-scale numerical simulations of 
stellar radiation zones ( |Garaud|2002[|Brun & Zahn|2Q06|| ). 



satisfies this type of surface boundary conditions. This will be 
the next step and it is out of the scope of the present paper. 



7. Conclusion 

In the context of improving stellar models by taking into account 
in the most consistent way as possible dynamical processes such 
as rotation and magnetic field, we examine possible magnetic 
equilibrium configurations to model initial fossil fields. 

We generalize the pioneer work by Prendergast (1956) in de- 
riving the barotropic magnetohydrostatic equilibrium states of 
realistic stellar interiors which are a first equilibrium family. 
These will then evolve due to other dynamical processes such 
as Ohmic diflfusion, diflferential rotation, meridional circulation, 
and turbulence. Relaxed minimum energy equilibrium configu- 
rations are then obtained for a given conserved mass and helic- 
ity that correspond to the Taylor's relaxation states in the self- 
gravitating non force-free case. These are then applied to the 
internal radiation zone of the young Sun and to the radiative in- 
terior of an Ap star on the ZAMS. Mixed poloidal and toroidal 
magnetic configurations, potentially stable in stellar radiation 
zones, are obtained. 

Now, we have thus to study the stability of these magnetic 
topologies. Moreover, the case of general baroclinic equilibrium 
states have to be studied (Paper II). 

These equilibrium configurations have then to be used as 
possible initial conditions for rotational transport processes in 
stellar radiative regions that will allow to study internal stellar 
MHD over secular time- scales. 
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6.3. Relaxed configurations and boundary conditions 

Let us now discuss the boundary conditions we choosed. Since 
equilibrium states are known to minimize the energy/helicity ra- 
tio, we follo w the proce dure established by Chandrasekhar & 
IPrendergast ([1958]) and Woltjer" (1959b) to constrain the arbi- 
trary functions of the magnetohydrostatic equilibrium. This pro- 
cedure, which minimizes the energy with respect to given in- 
variants of the system (and in particular the helicity), assumes 
the following boundary condition B - er = that leads to an az- 
imuthal current sheet due to the non-zero latitudinal field at the 
upper boundary (Bq (r^up^ ^) ^ 0)- This is a potential source of 
instability and in the case of our configuration we have to evalu- 
ate its eflTect on the stability (cf. |Bellan|2000l ). 

Next, in a stellar context, we have to allow open configura- 
tions as observed and thus to match the internal solution with an 
external multipolar one. It remains then to be seen whether the 
invariants are conserved, as they are in the confined case ( Dixon 
let al.|1989| ). 

Finally, independently from the chosen type of configuration 
(confined or matched with a multipolar external field), we have 
to search solutions that allow the continuity of the magnetic 
field and of the associated currents at the boundaries to cancel 
the possible induced instabilities. This leads to an ill-posed 
problem which must be solved in a subtle way (see Monaghan 
|1976[maitik ov 2009). In the present state of art, no solution has 
been derived that both minimizes the energy/helicity ratio and 
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